National  Bureau  of  Standards 


Uim 


. 


HITED  STATES 

kRTMENT  OF 

1MERCE 

3UCATION 


Lt»T0'cQ 


|>«ris  o*  * 


AUG 


18  1970 


NBS  TECHNICAL  NOTE  521 


U.S. 
ARTMENT 
OF 
OMMERCE 

National 
Bureau 
of 
Standards 


Present  State  of  the 
Classical  Theory  of 
Quantitative  Electron 
Probe  Microanalysis 


NATIONAL  BUREAU  OF  STANDARDS 

The  National  Bureau  of  Standards  '  was  established  by  an  act  of  Congress  March  3,  1901.  Today, 
in  addition  to  serving  as  the  Nation's  central  measurement  laboratory,  the  Bureau  is  a  principal 
focal  point  in  the  Federal  Government  for  assuring  maximum  applicatidn  of  the  physical  and 
engineering  sciences  to  the  advancement  of  technology  in  industry  and  commerce.  To  this  end 
the  Bureau  conducts  research  and  provides  central  national  services  in  four  broad  program 
areas.  These  are:  (1)  basic  measurements  and  standards,  (2)  materials  measurements  and 
standards,  (3)  technological  measurements  and  standards,  and  (4)  transfer  of  technology. 
The  Bureau  comprises  the  Institute  for  Basic  Standards,  the  Institute  for  Materials  Research,  the 
Institute  for  Applied  Technology,  the  Center  for  Radiation  Research,  the  Center  for  Computer 
Sciences  and  Technology,  and  the  Office  for  Information  Programs. 

THE  INSTITUTE  FOR  BASIC  STANDARDS  provides  the  central  basis  within  the  United 
States  of  a  complete  and  consistent  system  of  physical  measurement;  coordinates  that  system  with 
measurement  systems  of  other  nations;  and  furnishes  essential  services  leading  to  accurate  and 
uniform  physical  measurements  throughout  the  Nation's  scientific  community,  industry,  and  com- 
merce. The  Institute  consists  of  an  Office  of  Measurement  Services  and  the  following  technical 
divisions: 

Applied  Mathematics — Electricity — Metrology — Mechanics — Heat — Atomic  and  Molec- 
ular Physics — Radio  Physics  - — Radio  Engineering  - — Time  and  Frequency  - — Astro- 
physics L> — Cryogenics.2 
THE  INSTITUTE  FOR  MATERIALS  RESEARCH  conducts  materials  research  leading  to  im- 
proved methods  of  measurement  standards,  and  data  on  the  properties  of  well-characterized 
materials  needed  by  industry,  commerce,  educational  institutions,  and  Government;  develops, 
produces,  and  distributes  standard  reference  materials;  relates  the  physical  and  chemical  prop- 
erties of  materials  to  their  behavior  and  their  interaction  with  their  environments;  and  provides 
advisory  and  research  services  to  other  Government  agencies.  The  Institute  consists  of  an  Office 
of  Standard  Reference  Materials  and  the  following  divisions: 

Analytical  Chemistry — Polymers — Metallurgy — Inorganic  Materials — Physical  Chemistry. 
THE  INSTITUTE  FOR  APPLIED  TECHNOLOGY  provides  technical  services  to  promote 
the  use  of  available  technology  and  to  facilitate  technological  innovation  in  industry  and  Gov- 
ernment; cooperates  with  public  and  private  organizations  in  the  development  of  technological 
standards,  and  test  methodologies;  and  provides  advisory  and  research  services  for  Federal,  state, 
and  local  government  agencies.  The  Institute  consists  of  the  following  technical  divisions  and 
offices: 

Engineering  Standards — Weights  and  Measures  —  Invention  and  Innovation  —  Vehicle 
Systems    Research — Product   Evaluation — Building  Research — Instrument  Shops — Meas- 
urement Engineering — Electronic  Technology — Technical.  Analysis. 
THE  CENTER  FOR  RADIATION  RESEARCH  engages  in  research,  measurement,  and  ap- 
plication of  radiation  tb  the  solution  of  Bureau  mission  problems  and  the  problems  of  other  agen- 
cies and  institutions.  The  Center  consists  of  the  following  divisions: 

Reactor  Radiation — Linac  Radiation — Nuclear  Radiation — Applied  Radiation. 
THE  CENTER  FOR  COMPUTER  SCIENCES  AND  TECHNOLOGY  conducts  research  and 
provides  technical  services  designed  to  aid  Government  agencies  in  the  selection,  acquisition, 
and  effective  use  of  automatic  data  processing  equipment;  and  serves  as  the  principal  focus 
for  the  development  of  Federal  standards  for  automatic  data  processing  equipment,  techniques, 
and  computer  languages.  The  Center  consists  of  the  following  offices  and  divisions: 

Information  Processing  Standards — Computer  Information  —  Computer  Services  —  Sys- 
tems Development — Information  Processing  Technology. 
THE  OFFICE  FOR  INFORMATION  PROGRAMS  promotes  optimum  dissemination  and 
accessibility  of  scientific  information  generated  within  NBS  and  other  agencies  of  the  Federal 
Government;  promotes  the  development  of  the  National  Standard  Reference  Data  System  and  a 
system  of  information  analysis  centers  dealing  with  the  broader  aspects  of  the  National  Measure- 
ment System,  and  provides  appropriate  services  to  ensure  that  the  NBS  staff  has  optimum  ac- 
cessibility to  the  scientific  information  of  the  world.  The  Office  consists  of  the  following 
organizational  units: 

Office  of  Standard  Reference  Data — Clearinghouse  for  Federal  .Scientific  and  Technical 
Information  ' — Office  of  Technical  Information  and  Publications — Library — Office  of 
Public   Information — Office  of  International  Relations. 


'  Headquarters  and   Laboratories  at  Gaithersburg,   Maryland,  unless   otherwise   noted;   mailing  address  Washington.   D.C.   20234. 

-  Located  at  Boulder.  Colorado  80302. 

11  Located  at  5285  Port  Royal  Road.  Springfield,  Virginia  22151. 


UNITED  STATES  DEPARTMENT  OF  COMMERCE 
Maurice  H.  Stans,  Secretary 

NATIONAL  BUREAU  OF  STANDARDS    •     Lewis  M.  Branscomb,  Director 


NBS 


TECHNICAL  NOTE  521 

ISSUED  AUGUST  1970 

Nat.  Bur.  Stand.  (U.S.),  Tech.  Note  521,  17  pages  (Aug.  1970) 
CODEN:      NBTNA 


Present  State  of  the  Classical  Theory  of 
Quantitative  Electron  Probe  Microanalysis 


Kurt  F.  J.  Heinrich 

Spectrochemical  Analysis  Section 

Analytical  Chemistry  Division 

Institute  for  Basic  Standards 

National  Bureau  of  Standards 

Washington,  D.C.   20234 


NBS  Technical  Notes  are  designed  to  supplement  the 
Bureau's  regular  publications  program.  They  provide  a 
means  for  making  available  scientific  data  that  are  of 
transient  or  limited  interest.  Technical  Notes  may  be 
listed  or  referred  to  in  the  open  literature. 


For  sale  by  the  Superintendent  of  Documents,  U.S.  Government  Printing  Office,  Washington,  D.C.      20402 
(Order  by  SD  Catalog  No.    C  13.46:  521).      Price  30  cents 


PRESENT  STATE  OF  THE  CLASSICAL  THEORY  OF 
QUANTITATIVE  ELECTRON  PROBE  MICROANALYSIS* 

Kurt  F.  J.  Heinrich 

Analytical   Chemistry    Division 
Institute   for   Materials   Research 
National    Bureau   of  Standards 
Washington,    Do    C.     20234 

Although  the  foundations  for  a  procedure  of  data  reduction 
in  quantitative  electron  probe  analysis  have  not  been  changed 
for  several  years,  there  has  been  progress  in  the  choice  of 
expressions,  parameters,  and  constants.   A  brief  account  of 
recommended  expressions  and  procedures  is  given.   Reference  is 
made  to  the  Standard  Reference  Materials  of  Au-Ag  and  Au-Cu 
alloys  issued  for  electron  probe  microanalysis.   These  are 
especially  useful  for  investigating  the  application  of  correc- 
tion procedures. 

Key  words:   corrections;  electron  probe  microanalysis;  quanti- 
tative analysis;  x-ray  spectroscopy. 

INTRODUCTION 

In  his  doctoral  thesis,  published  in  1951  [1],  Castaing 
laid  the  groundwork  for  a  theory  of  quantitative  electron 
probe  microanalysis.   He  pointed  out  that  the  intensity  of 
characteristic  primary  radiation  generated  in  multi-element 
specimens  is  approximately  proportional  to  the  mass  fraction 
of  the  emitting  element.   He  therefore  proposed  that  the 
corrections  which  need  to  be  applied  to  measured  x-ray  inten- 
sity ratios  to  obtain  mass  fractions  ("concentrations")  were 
as  follows: 


* 

This  paper  is  to  appear  in  Mikrochim.  Acta,  in  German,  as  part 

of  the  proceedings  of  a  meeting  on  electron  probe  microanaly- 
sis and  metallography,  held  in  Vienna,  (Technische  Hochschule, 
Wien),  October  22-25,  1969. 


1.  An  absorption  correction,  which  takes  into  account 
the  absorption  of  x-rays  emerging  from  the  specimen,  and 

2.  A  fluorescence  correction  which  cancels  the  effects 
of  indirect  x-ray  generation,, 

Castaing  also  suggested  that  for  accurate  analysis  the 
proportionality  between  concentration  and  intensity  of  emis- 
sion may  have  to  be  considered  an  insufficient  approximation. 
Hence,  he  proposed  a  form  of  atomic  number  correction  to 
correlate  more  accurately  the  mass  fraction  with  the  measured 
intensity  of  x-ray  emission. 

This  "classical"  scheme  of  data  reduction  has  essentially 
remained  the  commonly  used  approach  to  quantitative  analysis. 
It  has  been  described  in  many  publications,  including  the 
National  Bureau  of  Standards  Special  Publication  298  (1968) 
[2],  which  contains  several  references  used  in  this  paper.   In 
the  recent  past  there  has  been  considerable  discussion  of  the 
choice  of  parameters  and  constants  to  be  employed  within  the 
framework  of  this  procedure.   The  problem  remains  if,  instead 
of  the  classical  procedure,  a  Monte-Carlo  method  or  a  transport 
equation  procedure  is  used  to  calculate  the  mass-fractions  of 
the  elements  being  measured. 

It  is  the  purpose  of  the  present  publication  to  give  an 
account  of  the  specific  choices  of  parameters  and  constants 
which  are  available  at  present  for  the  computation  of  correc- 
tions for  quantitative  electron  probe  microanalysis.   For  a 
general  description  of  the  correction  procedures,  as  well  as 
for  a  detailed  discussion  of  the  origin  of  these  expressions, 
the  reader  is  referred  to  the  cited  literature. 


ATOMIC  NUMBER  CORRECTION 
Thanks  to  the  work  of  many  investigators  [e.g.    2-6]  the 
processes  resulting  in  primary  x-ray  emission  are  now  quite 
well  understood.   The  number  of  ionizations  produced  by  a 
single  electron  coming  to  rest  within  the  target  can  be 
written  as  follows  ; 


"-"i/'i^  <» 


For   the    ionization   cross-section,    Q,    and   the   stopping  power, 
-dE/d(px),    we   can   substitute    in  the   above    equation  the   expres- 
sions  proposed   by   Bethe   et   at.    [7,8]: 

Q  =   b  •  n  e4  .  z£  .  Ml  .  E~2  (2) 

S.-i.i|S.W.NAT.|.It.    feV  <3) 

P     dx  Av     A    E        \J  / 

The  values  of  some  parameters  of  these  equations  are,  however, 
poorly  known „  Henoc  has  concluded  from  the  study  of  experi- 
mental results  on  the  ionization  cross  section  that  the  value 
b  =  0.76  ,can  be  considered  valid  for  all  levels  of  ionization 
[9].  The  value  of  the  coefficient  a  in  equation  (3)  was  fre- 
quently considered  to  be  equal  to  two.  However,  according  to 
Nelms  [10],  for  electrons  the  value  a  =  1.166  should  be  used. 

There  is  some  controversy  concerning  the  best  expression 
for  the  mean  ionization  potential,  J.  The  value  (J  =  11»5Z). 
used  frequently  in  the  past,  is  certainly  inaccurate.  We  use 
at  present  an  expression  obtained  by  Berger  and  Seltzer  [11], 
who  averaged  a  large  number  of  experimental  observations  by 
several  authors: 

J  =  9.76  Z  +  58.5  Z-°-19  (4) 


* 

A  list  of  symbols  used  is  appended. 


Another   expression   for   J   has   been   proposed  by  Duncumb   et   at. 
[12]: 


J/Z  =    14.0  [1-exp  (-0.1  Z)]    +   75.5/Zz/75-Z/(100  +  Z) 


(5) 


This  equation  has  been  derived  empirically  to  optimize  the 
results  of  a  large  number  of  electron  probe  analyses  of  speci- 
mens of  known  composition.   We  withold  judgment  on  the  best 
choice  of  expression  for  J,  pending  further  investigations 0 


In  order  to  obtain  an  algebraic  expression  for  ionization 
within  a  multi-element  target,  we  introduce  the  variables  M 
and  V,  defined  by  the  following  equations: 


=  2c3 

i  l  A. 


V 


exp 


M  i  Ci  A.  £n(!^ 


aE 
J 


(6) 


As  shown  by  Springer  [13] >  if  one  introduces  equations  (2), 
(3),  and  (6),  the  integration  of  equation  (1)  yields  the 
expression 


nA  =  I  JL  H   bIiu0-l-inY[EI(£n(U0V))-EI(£nV)]l 
2AM1       V  ' 

in  which  EI  is  the  exponential  integral: 


(7) 


EI  (x )  =  const.  +  £n  x  +  2 


8  =  1 


s-  s  ! 


(8) 


Expression  (7)  is  used  in  an  extensive  program  for  data  pro- 
cessing in  quantitative  electron  probe  microanalysis,  called 
COR,  recently  developed  by  Henoc  and  Heinrich  at  the  National 
Bureau  of  Standards  [14], 


A  further  correction  is  necessary  to  take  account  of  the 
effects  of  electron  backscattering.   A  fraction  of  the  imping- 
ing electrons  is  reemitted  with  energies  sufficient  to  excite 
x-rays.   This  reduces  the  generated  x-ray  intensity  by  a  fac- 
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tor  R  (backscatter  effect),  which  decreases  rapidly  with  in- 
creasing atomic  number.   The  average  number  of  ionizations  per 
impinging  electron  is  thus  n»R.   We  believe  that  the  most 
reliable  expression  for  the  factor  R  is  that  proposed  by 
Duncumb  [15] : 

R  =  1.0000  + 

+  (-.581+2. 162W-5.137W2+9.213W3-8.619W4+2.962W5)   x    10-2Z  + 

+  (-1.609-8. 298W+28.791W2-47.744W3+46.540W4-17.676W5)   x    10"4Z2   + 

+  (5.400+19.184W-75.733W2+120.050W3-110.700W4+41.792W5)   x    10-6Z3   +  (9) 

+  (-5.725-21.645W+88.128W2-136.060W3+117.750W4-42.445W5)   x    10~8Z4   + 

+  (2.095+8.947W-36.510W2+55.694W3-46.079W4+15.85lW5)   x    10~10Z5 

where   W   =    1/U      =   E   /E,  .      This    expression   agrees   well   with   the 

O      O   K! 

experimental  measurements  of  the  factor  R  by  Derian  [16], 
The  number  of  photons  of  the  x-ray  line  n,  generated  by  ion- 
ization of  the  electron  shell  m,  is  obtained  by  multiplying 
the  number  of  ionizations,  n.R,  with  the  fluorescence  yield, 

a)  ,  and  with  the  ratio  of  the  intensity  of  the  line  of  interest 
m'  J 

to  the  intensity  of  all  lines  originating  from  this  shell 

(weight  of  the  line,  p   ): 
&>  i    ^mn 

I'   =  nA  .  R.o)  .  p  (10) 


THE  ATTENUATION  OP  PRIMARY  RADIATION 

The  attenuation  of  primary  radiation  on  emergence  from 
the  specimen  can  be  calculated  using  an  equation  proposed  by 
Philibert  [17] 

i/f(x)  =  (i  +  £>(i  +  rr-r-  ■-£  >;  h  =  1-2  A.;      (11) 

The  factor  f(x)  is  the  ratio  between  attenuated  and  unattenu- 
ated  primary  radiation:   F(x)  =  In/I f •   Expressions  for  the 
coefficient  a,  which  depends  on  the  electron  energy,  were  pro- 
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posed  by  Duncumb  and  Shields  [18]  and  by  Heinrich  [19] •   We 
presently  use  the  following  equation: 


4.5  x  105 
E  !-65  -  E,  1-65 

o         k 


(12) 


We  have  not  evaluated  an  alternative  equation  for  the  absorp- 
tion correction,  by  Wittry  and  Andersen  [20]. 

The  attenuation  of  primary  radiation  is  large  when  x-rays 
of  low  photon  energy  are  measured.   It  is  important  to  note, 
however,  that  the  magnitude  of  this  effect,  and  hence,  the  un- 
certainty of  its  computation,  depends  upon  experimental  condi- 
tions under  our  control.   Primary  absorption  can  be  minimized 
by  working  at  low  electron  energies  and  using  high  x-ray 
emergence  angles;  and  these  conditions  are  desirable  when  soft 
x-rays  are  measured.   It  should  be  added  that  further  experi- 
ments on  the  depth  distribution  of  primary  x-rays  should  be 
encouraged.   The  uncertainties  in  the  absorption  correction 
have  been  analyzed  by  Yakowitz  and  Heinrich  [21]  who  suggested 
that  the  experimental  conditions  should  be  adjusted,  wherever 
possible,  to  maintain  the  value  of  f(x)  above  0.8. 


SECONDARY  RADIATION  EXCITED  BY  CHARACTERISTIC  LINES 


The  effects  of  secondary  radiation  excited  by  character- 
istic lines  are  usually  treated  according  to  the  equation 
proposed  by  Castaing  [1]: 


rf  .  h±  =  C 

'pa 


!pa  rA  -  *  wb  /*  <A,BK«)    A 


B  " 


XpA 


2  ^(AB,BK«)   B 


Ml  +  u) 


(1  +  v) 


(13) 


In  this  equation,  r   is  the  ratio  of  the  attenuated  fluorescent 

radiation  to  the  attenuated  primary  radiation  of  the  element 

which  is  being  determined.   I'  /I'   is  the  ratio  of  the  gener- 

pB   pA 


ated  (not  attenuated)  primary  emissions  of  the  exciting  and 
the  excited  lines  as  obtained  from  pure  elements „   Castaing 
proposed  to  approximate  this  ratio  by  the  ratio  of  the  absorp- 
tion edge  wavelengths  of  these  elements,  while  Reed  [22]  re- 
placed it  by  the  expression  [(UB  -  1)/(UA  -  D]  L67  in  which  U"A 

and  uV.  .are  the  corresponding  overvoltage  ratios. 
B 

This  procedure  is  amenable  to  several  improvements  which 
were  incorporated  in  the  program  COR.   In  this  program,  the 
ratio  I,-,/ In  is  fully  calculated,  using  equation  (10)  for 
primary  intensities,  just  as  in  the  atomic  number  correction. 
This  simplifies  the  calculation  of  fluorescence  of  K  lines  by 
L  lines  and  viae    versa,    as  well  as  the  calculation  of  fluores- 
ence  of  Ka  lines  by  K3  lines,  which  is  performed  separately. 
As  pointed  out  by  Criss  [23],  the  exponential  approximation  to 
primary  distribution  in  depth  of  x-ray  emission,  used  by 
Castaing  to  calculate  the  term  v  in  equation  (13) 9    can  be  re- 
placed by  a  more  accurate  model.   We  use  in  COR  the  model  con- 
tained in  Philibert's  calculation  for  the  function  f(x)  [17]. 
With  the  computer  facilities  presently  available  to  most 
analysts,  this  more  rigorous  approach  is  possible  without  an 
excessive  increase  in  the  cost  of  computation. 

SECONDARY  RADIATION  EXCITED  BY  THE  CONTINUUM 

Henoc  has  described  a  method  for  calculating  the  effects 
of  secondary  radiation  excited  by  the  continuum  [24].   In 
order  to  perform  the  formal  integration  of  the  expressions 
contained  in  his  procedure,  he  had  to  make  the  simplifying 
assumption  that  between  absorption  edges  the  mass  attenuation 
coefficients  vary  with  the  cube  of  wavelength.   However,  the 
wavelength  dependence  of  mass  attenuation  coefficients  is  al- 
ways less  than  cubic  [25].   This  has  been  taken  into  account 
in  the  modification  of  Henoc's  original  method  incorporated 
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in  the  program  COR.   The  formal  integration  is  replaced  by  a 
numerical  integration,  using  the  parameters  for  attenuation 
coefficients  published  in  reference  [25], 

The  main  practical  difficulty  in  performing  the  correc- 
tion for  fluorescence  excited  by  the  continuum  is  that  the 
position  of  the  absorption  edges  of  all  elements  contained  in 
specimen  and  standard  must  be  taken  into  account.   Each  wave- 
length range  between  adjacent  absorption  edges  must  be  calcu- 
lated separately.   Therefore,  the  procedure  is  lengthy,  and 
in  machine  computation  it  requires  extensive  data  input. 
Errors  in  the  input  occur  frequently  and  may  escape  detection, 
since  it  is  difficult  to  acquire  an  intuitive  feeling  for  the 
magnitude  of  this  correction. 

In  the  program  COR  the  relations  between  lines  and  absorp- 
tion edges  are  inspected  by  the  computer  which  selects  the 
wavelength  ranges  for  the  integration.   These  decisions  are 
made  and  the  computer  completes  the  calculation  without  active 
intervention  of  the  operator.   In  the  same  process,  the  com- 
puter also  determines  which  operations  must  be  performed  for 
calculating  secondary  excitation  by  characteristic  lines.   This 
enormously  simplifies  the  execution  of  the  fluorescence 
corrections,  and  renders  unnecessary  further  attempts  to 
simplify  them,  as  advocated  by  Springer  [26]. 

THE  ITERATIVE  PROCESS 

We  will  now  discuss  how  the  individual  corrections  are 
combined  in  calculating  the  composition  of  the  specimen. 

The  usual  form  employs  multiplicative  correction  factors: 

C*  =  k-  k*  .  k^  •  k*  (14) 


This  procedure  is  based  on  the  concept  that  ideally  the  method 
should  yield  a  linear  correction  curve.   The  multiplicative 
factors  correct  for  deviations  from  the  ideal  model.   However, 
the  values  of  the  correction  factors  depend  on  the  composition 
of  the  specimen,  which  is  unknown.   The  problem  is  tradition- 
ally resolved  by  using  an  estimate  of  composition  to  calculate 
the  correction  factors,  and  using  the  resulting  mass  fractions 
as  a  new  estimate  of  composition  and  iterating  until  converg- 
ence of  results  is  obtainedD 

As  stated  by  Criss  [23],  the  formulation  of  the  correc- 
tion by  multiplicative  factors  unnecessarily  complicates  the 
structure  of  the  correction  equations,  and  obscures  the  under- 
lying physical  events.   Moreover,  the  iteration  frequently 
fails  to  converge,  as  observed  by  Reed  [27].   For  these  reasons, 
the  program  COR  uses  a  different  iterative  concept,  first  pro- 
posed by  Criss  and  Birks  [28],  and  somewhat  modified  by  Henoc 
and  Heinrich. 

Each  iteration  contains  two  steps.   In  the  first  step,  we 
establish  an  estimate  of  the  specimen  composition  and  calculate 
the  intensity  ratios  to  be  expected  for  each  element,  according 
the  theory.   The  intensity  ratios  are  obtained  by  the  simple 
relation 

I*   I*  +  I*  +  SI* 

K  =-  =  (15) 

I    Ip  +  Ic  +  Sif 

The  correction  procedures  previously  discussed  are  used  to  cal- 
culate the  absolute  intensities  (in  photons  per  electron)  due 
to  direct  and  indirect  excitation.   As  a  first  estimate  of  the 
specimen  composition  we  assume  that  the  weight  fractions  of  the 
elements  are  proportional  to  the  experimentally  obtained  x-ray 
intensity  ratios.   This  estimate,  as  well  as  the  subsequent 
estimates  is  normalized  so  as  to  render  the  sum  of  mass  frac- 
tions equal  to  unity.   Thus,  in  each  iteration,  we  calculate 
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the  signal  intensities  to  be  expected  from  a  theoretically 
possible  composition. 

Since  equation  (15)  does  not  yield  an  explicit  term  for 
mass  fractions,  we  must  use  another  relation  to  obtain,  in  the 
second  step  of  each  iteration,  a  new  estimate  of  composition. 
We  use  for  this  purpose  a  hyperbolic  approximation  to  the 
analytical  calibration  curve  for  each  element: 

1-k     1-C  (1,. 
=  a  (16) 


Ziebold  and  Ogilvie  observed  that  the  analytical  curves  in 
electron  probe  microanalysis  are  indeed  approximately  hyper- 
bolic [29].   On  the  plot  of  intensity  ratio  versus   mass  frac- 
tion for  each  element,  the  hyperbola  can  be  defined  by  three 
points  (0,0),  (1,1),  and  (C,k),  where  C  and  k  are  the  estimate 
of  mass  fraction  used  in  the  iteration,  and  of  the  intensity 
ratio  calculated  for  it.   Our  next  estimate  for  the  mass 
fraction,  C  ,  corresponding  to  the  observed  intensity  ratio, 
k  ,  can  thus  be  calculated  by  the  equation 

c       km  ■  C(l-k)  (17) 

km(C-k)  \   k(l-C) 

Contrary  to  the  procedure  suggested  by  Criss,  we  do  not  norm- 
alize the  observed  intensity  ratios  k  in  equation  (16),  or  the 
final  results  after  convergence.   Thus,  if  the  presence  of  an 
element  has  been  overlooked,  or  if  there  are  experimental 
errors  or  failures  of  the  correction  procedures,  this  will 
become  apparent  since  the  calculated  mass  fractions  fail  to 
sum  to  unityc 

This  iteration  procedure  has  been  extensively  tested.   We 
have  never  observed  failure  to  converge.   In  fact,  the  converg- 
ence is  always  rapid,  as  is  to  be  expected  in  view  of  the  good 
fit  of  the  hyperbolic  model  observed  by  Ziebold  and  Ogilvie. 
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The  program  COR  can  also  compute  results  of  analysis  with 
use  of  multi-element  standards;  furthermore  it  can  determine 
one  element  by  difference  or  several  unmeasured  elements  on 
the  basis  of  stoichiometric  relations. 

STANDARD  REFERENCE  MATERIALS 

In  spite  of  the  progress  made  in  recent  years,  there  are 
still  uncertainties  in  several  aspects  of  the  correction  pro- 
cedures o   I  believe  that  the  availability  of  reliable  standard 
reference  materials  will  be  instrumental  in  pointing  out  the 
areas  in  which  further  investigation  is  most  needed.   For  this 
reason,  the  National  Bureau  of  Standards  has  devoted  consider- 
able effort  to  preparing  and  characterizing  materials  specially 
selected  for  electron  probe  microanalysis  [30].   Recently,  NBS 
certified  and  issued  an  alloy  of  nominal  composition  0.8  W, 
0.2  Mo,  a  series  of  binary  alloys  in  the  Au-Ag  system,  in  nom- 
inal steps  of  0.2  (20  weight  percent),  and  a  similar  series  of 
alloys  in  the  system  Au-Cu.   All  these  alloys  were  carefully 
analyzed  chemically  by  several  laboratories,  and  painstakingly 
investigated  for  homogeneity.   They  should  be  particularly  use- 
ful in  testing  the  atomic  number  correction,,   Extension  of  the 
studies  of  oxydic  systems  is  presently  in  progress. 
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List    of   Symbols 

A,B  atomic   weights    of   elements    A,    B 

C  mass    fraction    ("concentration") 

E  kinetic   energy    of  the   electron 

E  initial   kinetic   energy    of   the   electron,    operating 

voltage 

E,  excitation   potential    of   line   k 

k 

I  attenuated    (emergent)    intensity    of   radiation 

t 
I  intensity    of   generated   radiation,    before    attenuation 

J  mean   ionization   potential 

M  variable   used    for    stopping  power   calculation    (eq. 

6-7) 

Nn  Avogadro's  number 

Av 

Q  ionization   cross-section 

R  correction   factor   for   electron   backscattering 

S  stopping   power 

U  =   E/Ek 

U  =    E   /E" 

o  ok 

V         variable  used  for  stopping  power  calculation  (eq. 
6-7) 

W         =  1/U 

o 

Z  atomic   number 

a  constant    in  the   hyperbolic    approximation 

b  constant    in  the    calculation    of   the   ionization 

cross-section 

c  constant    in  the    calculation    of   a 

e  charge    of   electron 
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f(x)       correction  factor  for  absorption  of  primary  radiation 

2 
h  =    1.2(A/Z    )    in   Philibert's   equation.       For   multi- 

element   targets,    ,     _    y   r    .      .„    2 

i    i      i 

1 

EI  exponential   integral 

In  .       natural   logarithm 

nfl  number    of   ionizations   produced   by    one    electron 

which    comes   to    rest   within  the    specimen 

p  relative    frequency    of   the    line   n,    generated    on 

ionization   of   the    shell   m,    with   respect    to 
generation    of   all    lines   originating    from   shell   m 

r  absorption   edge    jump    of   the   absorption    coefficient 

r„  intensity    of   emerging   fluorescent    radiation, 

relative   to    emerging   primary    radiation    of   the 
same    line 

s  tally    in  the   exponential   integral 

k  relative    intensity    of   emergent    radiation 

k„  jk.  jk-p,      correction    factors    for    atomic   number,    absorption, 
and   fluorescence   effects 

u,v  auxiliary    variables    for    attenuation    of    fluorescent 

radiation 

a  constant    in   Bethe's   equation    for   energy    loss    of 

electrons 

6  emergence    angle    of  measured   x-rays 

y  x-ray   mass    attenuation    coefficient 

y(AB,AKa)    mass   attenuation    coefficient    of    specimen   containing 
elements   A   and   B,    for   Ka  radiation    of   element   A 

p  density    of   specimen 

a  coefficient    for   electron   energy    effect    in   Philibert's 

equation 

X  =    y    cosec    6 

oo  =    fluorescence   yield 
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Subscripts  and  superscripts: 

A,B       elements  A,  B 

i         tally  in  integration 

1         electron  shell  1 

*         composition  dependent 
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